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Abstract 

Medium modification of pion form factor has been evaluated in 
asymmetric nuclear matter (ANM). It is shown that both the shape 
and the pole position of the pion form factor in dense asymmetric 
nuclear matter is different from its vacuum counterpart with p-uj mix- 
ing. This is due to the density and asymmetry dependent p-oj mixing 
which could even dominate over its vacuum counterpart in matter. 
Results are presented for arbitrarily mixing angle. Effect of the in- 
medium pion factor on experimental observables e.g., invariant mass 
distribution of lepton pairs has been demonstrated. 

25.75.Dw: 21.65.+f: 14.40.-n 



1 Introduction 

The pion electromagnetic form factor, Ft^{Q'^) in vacuum shows that the 
physical p meson is not a pure isospin eigenstate and it can mix with the 
uj meson [1]. In vacuum this mixing amplitude can be determined by mea- 
suring Ft^IQ"^) which, although dominated by the p° pole, shows a kink near 
the u meson mass. Such mixing (after electromagnetic correction) implies 
that the charge symmetry is broken at the most fundamental level in strong 
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interaction through the small mass difference between up and down quarks 
in the QCD Lagrangian. Consequently, the physical p and oj mesons that 
we deal with are admixtures of the corresponding isospin eigenstates. At 
the hadronic level this mixing can be understood in terms of neutron-proton 
mass difference in effective models [2]. 

p-uj mixing has important and interesting consequences. It plays a cru- 
cial role in generating contributions to few body charge symmetry violating 
observables [3]. The p-uj mixing amplitude is determined from e^e~ — > tt'^tt" 
by measuring the pion form factor in the interference region [1]. With the 
extracted value of the mixing amplitude one is able to explain a number 
of observables namely, the non-Coulombic binding energy difference (Nolen 
Schiffer anomaly [4, 5]) of = 3 (mirror) nuclei, significant contributions 
to the np asymmetry at 183 MeV and non-negligible contributions to the 
difference of nn and pp scattering lengths [6] . 

The mixing of different isospin states will be modified in matter. Such 
medium effects have recently been investigated by several authors [7, 8, 9]. 
Unlike neutron-proton mass difference, which is responsible for p-uj mixing 
in free space, the mixing in matter can be induced if the neutron-proton 
densities are different. This happens even if the Hamiltonian preserves the 
isospin symmetry i. e., if M„ = Mp, akin to the 'spontaneous symmetry 
breaking' driven by the n p asymmetric ground state. It has been shown 
in Ref.[7] that the density dependent mixing is of similar magnitude as the 
usual vacuum mixing at normal nuclear matter density. Subsequently, Bro- 
niowski and Florkowski, showed that the mixing could be significantly large 
in Pb like nuclei [8]. These studies were limited to hadronic models. In con- 
trast, such density dependent p-u mixing was studied within the framework 
of QCD sum-rule in Ref. [10] with similar results and conclusions. The sum 
rule calculation has recently been improved and extended to study the den- 
sity dependent p-cu amplitude and its observable implications in heavy ion 
collisions [11]. 

In the present paper we revisit the problem within the framework of 
Walecka model [12]. We discuss the possible consequences in presence of 
a scalar mean field, an effective way to incorporate nucleon-nucleon inter- 
actions. Moreover, as in dense matter with large neutron-proton density 
asymmetry, the mixing angle could be quite high, we present results valid 
for arbitrary mixing angles. As an application, the possible modification of 
pion form factor in ANM and the Tr~^7r~ — > e~^e~ annihilation cross section at 
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various densities and asymmetries have also been discussed. This is relevant 
for the dilepton production in relativistic heavy ion collisions. Therefore, the 
present investigation has direct relevance for the study of compressed bary- 
onic matter (CBM) expected to be produced in heavy ion coUisions at GSI 
energies [13]. 

Apart from p — uj mixing in asymmetric matter p can also mix with the 
cr in the same way as a; — cr mixing [14] . We have discussed such a possibility 
in Appendix I where it is argued that for matter with small neutron-proton 
asymmetry such possibilities could be ignored. 

The paper is organized as follows. In the next section the formalism is set 
forth. In section 3 the relations between dilepton production cross section 
and form factor with mixing has been presented. Results are discussed in 
section 4 while section 5 is devoted to summary. Mathematical expressions 
relevant for the calculation of pion form factor with matter induced mixing 
and p and uj meson self energies in dense nuclear matter are relegated to the 
appendix. 



The Lagrangian describing p - uj meson-nucleon interaction is given by, 



Here M = {Mp+Mn)/2 where the subscript p (n) stands for proton (neutron). 
Clearly, due to the Pauli matrix ra, p meson couples to neutron with a 
negative sign in contrast to the cu meson while with proton both couple with 
the same sign. This gives rise to mixing that vanishes in the hmit, M„ = Mp 
which is illustrated in Fig.l. 

Here the mixing is generated by NN loop. The vacuum p-LJ mixing in 
this approach was first calculated in Ref.[2]. We calculate similar mixing in 
ANM which, besides n-p mass difference, would also depend both on baryon 
density (pb) and the asymmetry parameter a — {pn — Pp)/{Pn + Pp)- 



2 Formalism 
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Figure 1: p and uj nucleoli loop 



To determine the density and asymmetry dependent mixing amplitude we 
evaluate these loops at finite density which are characterized by the following 
density dependent polarization functions : 

= J^J'^'^ Tr[zr^zG(ii' + Q)ifliG{K)] (2) 



where, 



G{K) = Gf{K) + Gd{K) (3) 

G.(i^)Hi^.r + M*)[^,-^L__] (4) 

Gr>{K) = {K,r + M*)[—5{k, - El)e{kp - \k\)] (5) 

Here Gp denotes free nucleon propagator while Gd represents the density 
dependent part which forbids on mass shell nucleon propagation in matter 
due to Pauli blocking [12]. As the propagator now explicitly involves the 
Fermi momentum, contributions will be different for the neutron and proton 
loops if their densities (and hence Fermi momenta) are different. It is to be 
noted that the polarization function or the mixing amplitude will have two 
parts, one is like the polarization tensor with the free nucleon mass replaced 
by its effective mass called the free part representing nucleon-antinucleon 
excitations (Dirac sea) and the other one is density dependent part relevant 
for the scattering from the Fermi sphere. We discuss them separately. 
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2.1 Density dependent part 

The density dependent piece of the mixed polarization (p-cu) due to p-p or 
n-n excitations is generically given by 



(27r)4 



k\)x 



T^,{K-Q,K) 



{K - QY - M*2 

T^,{K,K + Q) " 



*2 



where T^y represents relevant traces as in Ref.[15]. 

The n^^(g) — 11^^ + IIJ^*^ functions in this case are as follows 



^ TT"^ Jo 



El Q^-A{K-Qy 



(6) 



(7) 



Ut{kF,M\qo,\q\ 



gpg^.{kp)M* 



8M 



)2Q^Q 



X 



d^k 



e;q^-a{k-q)^ 



(8) 



w here JC,, = {K, - ^Q^){K, - ^Q.), Q^. = {-g,. + ^) and E*, = 
\/k'^ + M*^. In the above expressions, for proton (neutron) loop we substitute 
M and kp with Mp (M„) and kp (kp) respectively. Moreover, at this point 
it might be recalled that for a vector meson moving in nuclear matter the 
longitudinal (L) transverse (T) polarization tensors are different because of 
)Cfj,i,, unlike the vacuum part which is proportional to Qf^^- The L and T 
modes are constructed as 11^ = — IIoo + 1^33 and 11^ = 1122 = I^ii, where the 
meson momentum Q = (go, 0, 0, \q\) (see Appendix for details). 

The mixing amplitude is characterized by H'lt which involves scattering 
from the neutron and proton Fermi spheres : 

ni> = n£>(fc^, m;, go, \q\) - n£;^(fc^, m^, go, \q\). (9) 
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The negative sign arises because of the in the p-NN interaction. 

The pure part of the polarization can be obtained by taking appropriate 
vertex factor like g^^ or Qp for both the vertices. Accordingly, the total is 
given by a sum over the neutron and proton loops instead of the difference. 

Kt = KTikF,M;,qo, \q\) + UZAk'^, qo, \q\) (10) 

and similarly p ^ uo gives results for the u meson polarization function. 
We take = 10.1, ^^^^ = Q and Qp = 2.63, Hp = 6.0 [16] in numerical 
computations. 



2.2 Free part 



The vacuum part will also give rise to mixing which is same as Ref . [2] with n 
and p mass M^^p replaced by the in-medium masses, M* ^. In Walecka model 
this is determined from the following self-consistent condition [12], 



M: 



n,p 



Mn,,-^yp + p'nl 



where p| {i = p, n) represent scalar densities given by 



Ml 

27r2 



E*kF - Mf In 



Mf 



The free part of the polarization tensor can be written as, 

n^'^ = (-^^'^ + Q''g7g')nF(g'). 

The mixing contributions to Hj? are given by : 



27r2 



dzz{l-z) 

Jo 



In 



Mf - z{l - z) 
M^-Q^z{l-z) 



(11) 



(12) 



(13) 



(p^n) 



(14) 



i-rpo; 
'-'-F,vt 
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m; 

Mr, 



f 

Jo 



dz\n 



m;^ - Q'' z(l - z) 
M'^-Q^z{l-z) 



{p — > n) 



(15) 
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The pure p {uo) meson self-energies for the vector-vector, vector-tensor 
and tensor-tensor parts are given by [16] 



pp{u!U!) 



27r2 



i=p,n 



+ 



M*L 



2Mi 

2 \ 2M; 



(Q2/f) + Mf/W) 



*2 7-(^)^ 



where 



ii) 



f 

Jo 



dz z(l 



rii) 



i: 



dz In 



z) In 

M*2 



M*' - Q' zjl - z) 
Mi-Q^z{l-z) 

-Q^z{l-zy 



(16) 



(17) 



(18) 



g2 z(i - z) 

It is to be noted that the free part 11^' vanishes in the hmit = Mp. 
We extract the real part of the vacuum mixing amplitude (with free nucleon 

mass to be 3447MeV2 at the omega pole and this is consistent with that 

of Ref. [17]. When in- medium nucleon masses are included 11 



put 



[m^j) is equal 

to —3716 MeV^ and —4675 MeV^ at po and 2po respectively in symmetric 
nuclear matter. 



3 Pion form factor and tt+tt — ^ e^e cross sec- 
tion 

The p-dominated (unmixed) pion electromagnetic form factor is given by: 

FAQ')-l-^-^yY, ?!■ ^ (19) 

in the first form of VMD [1]. 

In presence of mixing, the above expression for pion form factor should 
be replaced by (sec Appendix for detailed derivation) 

Ftv{Q^) = 1 - Qprnr ~ 3pnn tan € 

SpS^ LLp^ Qp 
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So 

X . 20 

Here tan2e = with Sp(^) = Q'^-m^^^^^ -^pp{^u) + imp(^)V ^.i^^). Here the 

couphng of the physical p and uj states to the photon is considered. In this 
form the uj — > 7r"'"7r~ decay is understood to proceed exactly like the p but 
modified by the mixing factor tan e and as appears in the denominator. 

In the small mixing limit the term quadratic in Tlp^ in Eq. 20 can be 
dropped and tane ~ e. Thus, to lowest order in the mixing parameter e, the 
above expression takes the form [1] 

sp 9p s^ 

The couphng constants used are g'j^^/An ~ 2.9, 5'p/47r ~ 2.0 [1] and 



9> 



UJ 



m^r{p e+e-) 



9p \ mpr{uj e+e-) 

= 3.5 ±0.18. (22) 

In matter, the value of gp^^T^ can acquire a density dependence due to the 
couphng with the A excitation as discussed in [18, 19, 20] in the leading 
density approximation. It was shown that the coupling could increase by 
~ 30 % near the p-pole at nuclear saturation density for finite width of the 
A. In the zero width approximation it could be as large as double near the 
p-pole at saturation density. However, as far as our results are concerned 
there will be no appreciable change in the dilepton invariant mass spectrum. 
The problem of density dependence of the vertex function is an interesting 
problem by itself and a systematic approach would be to evaluate the full 
p-spectral function in the presence of mixing. For the present purpose we 
ignore such effects in order to bring out the effect of mixed polarization tensor 
into clearer focus. 

The cross section for dilepton production from pion annihilation is inti- 
mately connected to the density dependent pion form factor which is given 

by, 



4 (T / 47T?2 

a(qo, \q\,p,, a) = ^^-^f - -^WMo, \q\)?- (23) 
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It is to be noted that unlike in vacuum, the pion form factor in nuclear 
medium depends both on qq and \q\. The dilepton emission rate in terms of 
the above cross section is given by 

m ^te.M.f.") M^Ti..(M/T)(l-ig) (24) 



dM 

where Ki is the modified Bessel function of second kind. 

It should be mentioned that in deriving the expression for in-medium 
-^7r(Q^) (Eq.20), the pions have been considered to be on shell. In relatrivistic 
heavy ion collisions the main contribution to dilepton production comes from 
pion annihilation where the pions are generally assumed to be on their mass 
shells. However, in a medium the pion dispersion relation is different from 
that of vacuum due its interaction. It would be interesting to extend the 
above calculation to incorporate this feature. This however, goes beyond the 
scope of the present work. 



4 Results 

In Fig. 2, the in-medium pion form factor for |g| = 200 MeV together with 
its vacuum counterpart is shown. In matter the pole position shifts towards 
lower invariant mass indicating the decrease of p and uj meson masses in 
matter. It might be mentioned that in medium, the mass modification is 
caused by two different mechanisms, viz. the scattering from Fermi sphere 
and excitation of the Dirac vacuum. While the former gives rise to an increase 
of their masses, the latter dominates resulting in an overall reduction. Near 
the pole, the mixing amplitude increases by large factors depending upon the 
value of the asymmetry parameter a. Clearly the density and asymmetry 
parameter dependent mixing is much larger than the mixing due to n-p mass 
difference alone. In Rcf. [11] no significant effect of matter induced mixing 
on was observed. However, in the present model this is found to be 
substantial which can be attributed to the tensor interaction. Moreover, we 
clearly show that the small angle approximation are not a good assumption 
in calculating pion form factor even for asymmetry relevant for heavy nuclei 
hkc Pb. 

Fig. 3 shows results for |g| = 600 MeV. Though the qualitative features 
remain similar, one can see that the mixing amplitude depends strongly on 
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Figure 2: Pion form factor as a function of invariant mass ( M — y/Q^) 
with mean field including both vacuum and Fermi sea contributions. The 
dotted and dot dashed line represent pion form factor in vacuum without 
and with mixing respectively. The dashed and solid line depicts the same 
for a = 0.3, and p = po in small angle approximation and for the arbitrary 
mixing angle respectively. The dot-dash-dashed line correspond to the case 
for arbitrary mixing angle at p = 2po- "Full (lowest order)" corresponds to 
Eq.(20) (Eq.21) 



the three momentum of the moving meson. In fact, the mixing is stronger 
for a slower vector meson. 



In Fig. 4, we show the strong asymmetry dependence of the mixing which 
shifts the pole of the pion form factor towards lower invariant mass region. 
Moreover, we also find that the results are quite sensitive to the asymmetry 
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Figure 3: Same as Fig.2 with q = 600 MeV. 



parameter {a). 

As an application of the density and asymmetry parameter dependent 
pion form factor, we calculate the pion annihilation cross section in nuclear 
matter. 



10 





' 1 


1 




1 ' ■ 




(Full mixing p^) 






0=0.7 




— Lowe si order (p ) 




/■ 


q=600 MeV 




— Vacuum mixing 




//" 






— No mixing 




f 






-■■ 




t 








/ 1 
/ 1 


I 

I 
i 
i 
i 


V, 

V 

\ 
\ 
\ 






/ \i ' 




\ 






/ \' ' 
^ 1''* 




\ 

\ 














1 









0.2 0,4 0.6 0.8 



M(GeV) 

Figure 4: Same as Fig. 3 for a = 0.7. 
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The dilepton production cross section is directly proportional to the 
square of the pion form factor and bears similar qualitative features. This 
is shown in Fig. 5. For completeness, we also present results for the dilepton 
production rates for various combinations of density and asymmetry param- 
eter at T = 100 MeV in Fig. 6. Evidently the medium modified pion form 
factor leads to enhanced production of dileptons in the low invariant mass 
region. 
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M (GeV) 



Figure 6: Dilepton production rate with and without mixing at two times 
the normal nuclear matter density. 



5 Summary 

In the present paper we have calculated pion form factor in asymmetric 
nuclear matter within the framework of relativistic mean field theory valid 
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for arbitrarily large mixing angle. It is known that the vacuum pion form 
factor shows a shoulder like behaviour near the p meson pole due to its ad- 
mixture with the isoscalar uj meson indicating isospin symmetry violation. 
In this paper, we have discussed the possible enhancement of such mixing 
in asymmetric nuclear matter where unlike the vacuum case, this symme- 
try breaking is driven by the ground state due to the difference of neutron 
and proton densities. Automatically this would modify the pion form factor 
as we have shown. Another interesting aspect is the shift of pole position. 
This is related to the presence of scalar mean field leading to the reduced 
nucleon mass in matter. In-medium pion form factor naturally influences 
the pion annililation cross section which proceeds through an intermediate 
p meson. Relativistic heavy ion collisions at GSI energies offers the unique 
opportunity to probe the in-medium pion form factor through dilepton mea- 
surements. Hence we have also calculated dilepton production rate due to 
pion annihilation in matter. In the low mass region the dilepton yield is 
enhanced due to dropping of vector meson mass as well as p-u mixing. 

We conclude with the comment that density and asymmetry dependent 
p — UJ mixing is an interesting problem related to charge symmetry violation 
in matter and it could be worthwhile to extend the formalism developed in 
the present work to other approaches such as chiral perturbation theory and 
NJL model which have explicit chiral symmetry. 



Appendix - I 

Expressions for mixing angles and pion form 

factor 

In this section we derive the expression for mixing angle and pion form factor 
for arbitrary mixing. In p — uj — a model because of the possibility of scalar- 
vector mixing the mixed polarization tensor has the following form in ANM, 











-1 












11(70; 


n, 



pa \ 

.a ■ (25) 



The propagator in the presence of mixining is the solution of 

L'-i = L>o-^-n, (26) 
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where Dq = d\Sig{dp, d^, d^) with rfj = 1/(Q^ — mf + ie),i = p, u, a. 
The mixed propagator can be obtained by matrix inversion: 



D 



I 





^pa^auj Sfj^ipuj 


^puj ^Lua ^uj'^pa 
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p^uja ^a^ujp 
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A 

^pa^LO p SpTlua 
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A 
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A 



(27) 



where A 



The polarization functions in D correspond to the mixing of the pure 
isospin cigcnstates, p/,c<j/ and cx/. This matrix can be expressed in terms of 
the propagators of the physical p, uj and a fields which is an admixture of the 
pure states. This is achieved by appropriate field rotation involving various 
mixing angles in isospin space : 







(pi 






UJj 









(28) 



where 5* is 3 x 3 mixing matrix involving density dependent mixing angles. 
Let us first compare the polarization tensors Ui^^ and Up^ which have the 
same three momentum dependence. Now since Ui^^ is a function of sum 

of the proton and neutron densities and Up^ is that of the difference, thus 
^uju » hlpcr for small values of the asymmetry parameter a. Taking Up^^ a; 
the propagator matrix takes the form: 

/ ^ 



D = 



\ 



Al 



Al 



Al 

SnSuj lln 



(29) 



Al Al 

where Ai = SpS^s^ - s^p^U^p - SpII^^n^^ 

Now it is also known that the scalar-vector mixing vanishes in the limit 
of vanishing 3-momentum and is very small for low momentum. However, 
Upuj remains non-zero even for ^ = 0. It is also to be noted that medium 
dependent mixing effects are dominant only in the low momentum region 
which allows us to put U^^^ ps; in this limit [14]. Under these assumptions, 
D reduces to block diagonal form in which a decouples from p and cu so that 
we have. 
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D 







\ 




(30) 



Consequently, S takes the form 



cose 


sin e 





— sin e 


cos e 











1 



(31) 



We henceforth work in the 2x2 subspace of p and a;, i. e. we use : 



cos e sm e 
— sin e cos e 



(32) 



The mixing angle can be deduced by diagonalising the vector meson propa- 
gator [1]: 





SujSp ^pcj 


spSuj n^j^ 




Sp 




SpSu — IIp^j 



s 





Sp 



from which we obtain 



tan 2e 



spSuj 



2n 



/!X<J 



5p "5a 



In the small mixing angle limit the above equation reduces to 



e = 



(33) 



(34) 



(35) 



For the derivation of pion form factor for arbitrary mixing angle we need to 
evaluate the matrix element for the process 7 — > tttt, which can be written 

as 



•' '■7— >7r7r 



. •' 'p/— >7r7r 



o)ss- 



SpSu) Hp^^ 



(36) 
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where it is assumed that T„ 



0. The physical amphtudes can be ob- 



tained from the above equation which follows : 



M 
M 



7 — >uj 



cos eM.y^pj — sineA^^_^i^, 



sin tM. 



(37) 




Figure 7: Contribution to p-uj mixing to pion form factor. 



Thus the amplitude for 7 — tttt can be written as 



•' '7— >7r7r 



+M^' 



xtaneA^ 



7— >a; 



(38) 



With these it is very easy to read out pion electromagnetic form factor from 
Fig. (7), 



1 - fi'pT. 



X- 



^ 9pmT tan e 



which in the small angle limit reduces to Eq.(20). 



(39) 
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Appendix - II 

In this appendix the mathematical expressions for polarization functions are 
presented. 
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